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It is well known that, for bright gamma-ray pulsars with high statistics above a few GeV, the phase 
averaged spectral energy distribution (SED) is harder than a simple exponential cutoff above the 
break. We perform phase-resolved spectral analyses of bright gamma-ray pulsars and demonstrate 
that, even over narrow phase ranges, the SEDs of gamma-ray pulsars above the break energy are 
harder than a simple exponential cutoff. We argue within a radiation-reaction limited curvature 
framework that this is indicative of non-stationary emission or emission from multiple zones. 
Eurther, we address a common problem faced when fitting hard spectral tails with a power-law 
times a sub-exponential function. Namely, that the sub-exponent parameter does not describe any 
parameters of physical models of pulsar emission. We introduce a simple analytical fit function 
to solve this problem. 
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1. Introduction 

The most commonly accepted mechanism to explain the emission from gamma-ray pulsars is 
curvature radiation, which is radiation from a relativistic charged particle moving along a curved 
magnetic field line. While the particular details of specific magnetospheric models may vary con¬ 
siderably, it is generally agreed that if the emission originates in the magnetosphere, particles are 
accelerated in regions where the plasma density is lower than the Goldreich-Julian density. These 
regions, referred to as gaps or emission zones, can support large electric fields which energise 
particles. These energetic particles then lose energy via curvature radiation and other processes 
[6, 9, 13, 14]. Magnetospheric models based on curvature radiation can generally reproduce the 
gamma-ray emission from pulsars observed by Fermi-1 .AT. Other models, which place the site of 
emission in the current sheet beyond the light cylinder are also reasonably successful at reproducing 
the observational data [11]. 

Within the curvature radiation framework, the spectrum of the emitted gamma rays follows 
a power-law with a spectral index related to the energy distribution of the emitting particles [8]. 
Beyond some critical energy, the acceleration gains of an electron are equaled by the radiative 
losses (the so-called radiation-reaction limit), yielding a gamma-ray spectrum that is exponentially 
suppressed [10]. Thus, for a single stable electron population emitting curvature radiation, the 
gamma-ray spectrum takes the form of a power-law times an exponential cutoff (PLEC). 

Within an outer magnetospheric emission framework, gamma-ray emission beams are ex¬ 
pected to by wide due to the large radius of curvature of the magnetic field along which the particles 
move, and due the fact electrons and positions are accelerated in opposite directions, leading to anti¬ 
parallel emission beams. Relativistic effects due to the large co-rotation velocity and the proximity 
to the neutron star distort the emission beams. These effects, coupled with time-of-fight differences 
between individual emission sites combine to create phases where emission beams converge (caus¬ 
tics) or diverge. In this way, the geometric framework of outer magnetospheric models directly 
produce the structure (peaks, bridges, off-regions, etc) seen in pulsar emission profiles. 

Given these considerations, in an outer gap curvature emission model, there is likely no 1-to-l 
mapping between the phase of an observed photon and the emission region within the magneto¬ 
sphere where the photon originated. If this is the case, then one should expect to see evidence for 
the emission from multiple zones within a phase-resolved pulsar spectrum. In this work we show, 
through phase-resolved spectral analyses of bright gamma-ray pulsars, that the spectral shape in 
a given narrow phase range is more compatible with a power-law times a sub-exponential cutoff 
(PLSEC) than a PEEC. We argue that this PESEC shape results from a sum of PEEC spectra and 
that this is evidence for multi-zone emission. Eurther, we introduce a simple functional form which 
is informed by these arguments and show that it can replace the PESEC form in pulsar spectral 
fitting. 

The remainder of this manuscript is structured in the following way. In Section 2, we explain 
how we performed our phase averaged and phase resolved spectral analyses. In Section 3, we 
discuss the results from our analysis and show that over narrow phase ranges, the emission from 
gamma-ray pulsars is harder than a simple PEEC and better described by a PESEC. In Section 
4, we explain how this is evidence for multizone emission or nonstationary emission zones and 
introduce a fit function intended to replace the PESEC in pulsar spectral fitting. 
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2. Fermi-LAT Data and Analysis Steps 

The analysis presented uses 5 years of Pass 7 Reprocessed photon data recorded by the Fermi- 
LAT. yonrce-class events falling within a 20 degree region-of-interest (ROI) are selected around 
each pulsar. A zenith angle cut of 100 degees is imposed on all events in order to minimize con¬ 
tamination from gamma-rays created in the Earth’s limb. Phase values are assigned to each photon 
using the Fermi TEMP02 plugin and pulsar timing models from Dr. Matthew Kerr^ [12]. All spec¬ 
tral analyses (phase-averaged, phase-resolved, and energy-band fitting) are performed by applying 
three iterations of the maximum likelihood fit. In each step the fit-tolerance condition is increased 
and sources with a low test statistic at are removed from the fit model. To ensure that final hi model 
mafched fhe dafa, residual maps are generated by subtracting a skymap derived from the fit model 
from a counts map derived from the data. These maps are inspected to ensure the scatter in the 
residuals is normal and that no large scale spacial structures or anomalies are present. 

Phase averaged analyses are performed first, starting with input fit models derived from the 
3FGE catalog^ [3]. The spectral and normalization parameters of the pulsar under study are free to 
float, as are a the normalization parameters of all sources within 4 degrees of the source. All other 
parameters are fixed fo fhe 3FGE value, bar fhe normalization of fhe galacfic diffuse model which 
is also free fo floaf. The besf-fif model from fhe phase-averaged fif is used as fhe inpuf model for 
fhe phase-resolved fit. For this analysis all parameters, bar the those associated with the pulsar, are 
fixed fo fhe value derived from fhe phase-averaged fil^. In each phase range, fhe pulsar specfrum is 
models as bofh a PEEC function 

( 2 . 1 ) 

(wifh fhe sub/super-exponenf b = \) and a PESEC function (wifh b free). The phase gafes for fhe 
phase-resolved analyses are chosen fo equalise fhe number of counfs from fhe pulsar in each phase 
bin. 

To make a specfral energy disfribufion (SED) from fhe phase averaged analysis, our dafa is 
parfifioned info 12 evenly spaced logarifhmic bins in energy wifh 4 bins per decade. On each of 
fhese energy bins, fhe fhree sfage analysis process is applied, modeling each source wifhin 4 degrees 
of each pulsar as a power law wifh specfral index fixed fo 2. Using fhe oufpuf of fhis analysis in 
each energy bin, fhe energy flux for fhe pulsar is calculated. 

3. Phase Resolved Spectral Shape of Each Pulsar 

Figure 1 plofs fhe phase-resolved specfral parameters of fhe Geminga and Vela pulsars under 
bofh fhe PEEC and PESEC model assumptions. In fhe PEEC fifs, fhe values for fhe specfral index 

'https://confluence.slac.stanford.edu/display/GLAMCOG/LAT+Gamma-ray+Pulsar+Timing+Models 

^An additional step is performed in the analysis of the Vela pulsar since the extended source, Vela-X (3FGL J0833.1- 
451 le), is spatially coincident with the Vela pulsar. The parameters of this source are derived by gating the data in the 
phase-range 0.8-1.0 where there is negligible emission from the Vela pulsar, and performing a likelihood fit. The model 
parameters of Vela-X are fixed in all subsequent analyses to the value returned from this pulsar-gated likelihood analysis, 
with the necessary scaling applied to the prefactor parameter. 

^The prefactor parameters were scaled to the size of the phase-gate in question in each fit. 


3 




On the Spectral Shape of Gamma-ray Pulsars Above the Break Energy 


Christopher Bochenek 




Figure 1: The phase-resolved spectral parameters of the Geminga pulsar (top) and Vela pulsar 
(bottom). The shaded histogram in each figure shows the phasogram of the corresponding pulsar. 
Note that for Vela, the range on the x-axis is restricted to be between 0 and 0.7. The red circles 
correspond to the spectral index in the PLEC fit, the blue squares to the PLEC break energy, and 
the red triangles to the PESEC subexponent. The red horizontal line indicates an index value of 
unity. 


and break energy in each phase bin are consistent with those obtained from the Eermi collabora¬ 
tion’s analysis of both pulsars [1, 2]. Erom the phase-resolved PESEC fits to Geminga, we note 
that the returned b parameter consistently takes a values less than one, and is only compatible with 
unity in 3 out of the 30 phase bins within the 1-sigma error bar. In the Vela analysis, the b parameter 
was found to be incompatible with unity in all tested phase gates. We also note that the test statistic 
values for each pulsar were higher when modeling with a PESEC than when modeling with a PEEC 
in each phase range. These results show that the phase-resolved spectra are better described by a 
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PLSEC than by a PLEC. Eor Geminga, bar the 3 phase regions where b is compatible with unity, 
we rule out b= \ confidence levels between 1.3a and 12.5a, with an average deviation of 3.7a. 
Eor Vela, we rule out ft = 1 at confidence levels befween 6a and 29a, wifh an average deviafion of 
10a. 


4. Discussion 


Abdo ef al. [2] and Celik & Johnson [5] have shown fhaf sub-exponenfially suppressed specfra 
are easily produced by summing fogefher several pure exponential shapes wifh a range of cufoff 
energies. This, fhey argue, explains fhe sub-exponenfial cufoff seen in fhe phase-averaged specfra of 
some pulsars, since if is fhe summation of fhe phase-resolved specfra which show significanf phase- 
dependenf variafion in fhe cufoff energy (see Eigure 1). We have shown here fhaf, for brighf gamma- 
ray pulsars, fhe individual phase-resolved specfra are fhemselves nof exponentially suppressed, buf 
are heifer described by sub-exponenfial breaks. 

Wifhin fhe confexf of outer magnefospheric curvafure emission models, we can infeipref fhe 
appearance of sub-exponenfial cufoffs in fhe phase-resolved specfra as evidence fhaf fhe emission 
observed af a given phase originales from several differenl particle accelerafion zones, each wifh a 
differenl radiafion-reacfion energy limif. Given fhe wide emission beams and fhe various relafivis- 
fic and lime-of-fighl effecls expecfed in oufer-gap curvafure emission, fhe convergence of differenl 
beams al fhe same phase is expecfed and should resull in fhe lype of sub-exponenfial speclral sig- 
nalure we have observed. Eeung el al. [7], have argued fhaf sub-exponenfial cufoff shapes appear 
in pulsars because fhe accelerating volfage in a given gap is unslable and Ihus fhe lime-averaged 
emission from a single emitting zone is a summation over various acceleration states. Either in¬ 
terpretation is compatible with the data. We note that for Geminga, the phases with the roundest 
cutoff (where the parameter b takes its smallest value) occurs in the leading edge of the first peak 
and the trailing edge of the second peak while the bridge region, in general, has the sharpest cutoff 
(b closest to unity). A similar behavior is seen is Vela. If the emission profiles are inferprefed 
wifhin an outer gap formalism, fhe phase-dependenf behavior of fhe b paramefer may help shed 
lighf on fhe emission geomefry. 

The clear appearance of sub-exponenfial cufoffs in phase-averaged and, now, phase-resolved 
pulsar specfra indicafe fhaf fhe PESEC funcfional form is increasingly required for pulsar speclral 
fitting. However, fhe b paramefer in fhe PESEC function has no physical analog in Iheories high- 
energy emission and if is highly degenerafe wifh fhe cufoff paramefer {Ec). These facls mean fhaf 
inleipreling fhe PESEC til paramelers wifhin lypical emission models is impractical. As discussed, 
fhe rounded sub-exponenfial shape seen in pulsars is attributed fo fhe sum over several purely 
exponenlial cutoffs wifh a range of differenl break energy. Wifh Ihis in mind, we inlroduce fhe 
following ad-hoc formula: 


dE~N^^\ J 


E' = lo(“+x/^) 


(4.1) 


which we call SUMPEEC. This formula is a summation of N PEECs, wifh fhe break energy of each 
ranging from fhe value 10“ MeV lo 10“+^ MeV in N logarilhmically spaced sleps. The values of 
A, r, a and j3 are free fo floal. Eq is fhe decorrelafion energy. The denominator of fhe term in 
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(a) Geminga (b) Vela 


Figure 2: Phase averaged spectral energy distributions of Geminga (left) and Vela (right) that shows 
the best fit lines from a PLEC, PLSEC, and SUMPEEC. 


the summation has effect of normalizing the area under each PEEC with the Eiow parameter fixed 
fo fhe value of fhe lower limif of fhe filled energy range. In fhe case, V = 1 and jS = 0, Ibis form 
reduces fo fhe simple PEEC. Using Ihis function lo til fhe phase-averaged specfrum of each pulsar, 
we found lhal increasing fhe value of N asympfolically improved fhe fil, wilh fhe besl-fil values 
remaining unchanged once N > 5. To produce Eigure 2, we fixed N lo fhe value 20. Inferpreled 
wilhin an ouler gap curvalure emission framework, fhe values of 10“ MeV and 10“+^ MeV specify 
fhe range of culoff energies presenl in fhe magnelosphere. Eor Geminga fhe cutoff energies span 
fhe range 1.22 ±0.11 to 5.1 ±0.2 GeV. Eor Vela, fhe cutoff energies span fhe range 1.35 ±0.13 fo 
9.8±0.5GeV. 

Eigure 2 clearly shows lhal fhe SUMPEEC funclion reproduces a shape very similar lo fhe 
PESEC and lhal if can be used lo fil pulsar speclra. By conslruclion, Ihe parameters of Ihe 
SUMPEEC funclion are easily interpreted wilh outer magnelospheric emission models and, Ihus, 
we argue if can replace Ihe PESEC funclion in pulsar speclral filling. 
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